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ABSTRACT 

Observations of minute-scale flares in TeV Blazars place constraints on particle accel- 
eration mechanisms in those objects. The implications for a variety of radiation mech- 
anisms have been addressed in the literature; in this paper we compare four different 
acceleration mechanisms: diffusive shock acceleration, second-order Fermi, shear accel- 
eration and the converter mechanism. When the acceleration timescales and radiative 
losses are taken into account, we can exclude shear acceleration and the neutron- 
based converted mechanism as possible acceleration processes in these systems. The 
first-order Fermi process and the converter mechanism working via SSC photons are 
still practically instantaneous, however, provided sufficient turbulence is generated on 
the timescale of seconds. We propose stochastic acceleration as a promising candidate 
for the energy-dependent time delays in recent gamma-ray fiares of Markarian 501. 
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1 INTRODUCTION 

Rapid TeV variability has recently been observed by the 
HESS and MAGIC teles copes in two blazars; P KS 2155-304 
jAharonian et all 120071 ) and Markarian 501 (|Albert et all 
I2OO7I ). The timescales of a few minutes are shorter, by at 
least an order of magnitude, than the light crossing time of 
the central black hole with a mass of order M = 10^ M©. 
This suggests that the variability is associated with small 
regions of the highly relativistic jet rather than the cen- 
tral region, alth ough the latter possibi lity cannot be com- 
pletely excluded. iBegelman et al.l l|2008l ) have examined the 
constraints that these observations place on the size and lo- 
cation of the emitting region. With an observed variability 
timescale ivar and a jet Lorentz factor T, the flare, if it moves 
out with the flow, occurs at a distan ce greater than rfvarT^. 
For r 50 (|Begelman et al~ and minute timescale 

variability this places the flaring region at a distance in ex- 
cess of one hundred Schwarzschild radii (rs) from the central 
black hole. The relativistic particles that are ultimately re- 
sponsible for the emission are then required to have been 
ejected from the central region along with the jet, and sub- 
sequently survived out to 100 rg where they then radiate 
away their energy quickly, or alternatively the particles are 
accelerated within the jet itself, close to the emission region. 
In this paper we examine the latter possibility and compare 
the predicted acceleration timescales for four acceleration 
mechanisms with the requirements for producing a minute 
timescale flare. 

There are three fundamental questions that we address 
in this paper. Firstly, we want to see whether particle ac- 



celeration could be considered instantaneous on the flaring 
timescale; is it realistic to assume a fully-developed power- 
law injection spectrum at the onset of flaring? Related to 
this question is whether we can exclude particular acceler- 
ation mechanisms on the grounds that they are too slow. 
Secondly, we are interested to see if we can find a set of 
parameters for which only cooling would be present. If we 
assume a given input spectrum for particles in a blob trav- 
elling with the jet, can we expect further acceleration to be 
always present? Alternatively, is it more likely that within 
minute to hour timescales particles would only lose their 
energy instead of being able to continue accelerating on the 
same timescales? Finally, can any of the mechanisms work 
on timescales comparable to the lags between lower- and 
higher-energy radiation observed in some objects, and what 
mechanisms could explain the acceleration of particles on 
the timescale of minutes? 

In Section 2 we discuss four possible acceleration mech- 
anisms and their associated timescales. These are (i) first- 
order Fermi acceleration, (ii) second-order Fermi, (iii) accel- 
eration in a shear fiow and (iv) the converter mechanism. 
The role played by radiation mechanism and loss timescales 
are then compared with the acceleration properties in Sec- 
tion 3. 



2 ACCELERATION TIMESCALES 

When particles are scattered by magnetic fluctuations, they 
gain energy whenever two subsequent scattering centres are 
moving toward each other, leading to a "head-on" colli- 
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sion. Suitable conditions are provided around a shock wave, 
where a relativistic particle crossing the shock always sees 
the plasma - and thus the scattering centres - of the flow on 
the other side of the shock, approaching. Similarly, relative 
motion of the scattering centres (Alfven waves, for exam- 
ple) inside a constant-velocity plasma, moving in different 
directions, are more often seen as approaching than reced- 
ing, leading to average gain in energy. The average energy 
gain during a crossing cycle, and the duration of such a cy- 
cle, determines the energy gain rate, < A'y/At >, and the 
acceleration timescale r = 7/ {d'y/dt). 

The timescale determines how long it takes for a par- 
ticle with Lorentz factor 7 to gain or lose energy, allowing 
easy comparison between the efficiencies of different gain 
and loss mechanisms. For easily comparable and general re- 
sults we assume the acceleration to take place in a region 
with a characteristic, comoving spatial scale R, with plasma 
flow having Lorentz factor F. Particles are assumed to be 
scattered by irregularities in large-scale magnetic field (with 
strength B) with the average free path between scattering, 
A, being equal to the particle gyroradius r-g. 

2.1 First-order Fermi acceleration 

The first-order process takes place at the shock front, where 
the accelerating particles gain energy by crossing and re- 
crossing the shock. An average cycle increases the particle 
energy by a factor of F^ for the first cycle, and by a factor of 
~ 2 thereafter. The duration of such a cycle, as well as the 
probability for a particle to be injected into one, depends 
heavily on the details of the scattering of the particles in 
the turbulent plasma and the geometry of the shock. In the 
Bohm limit where the particle's mean free path is equal to 
its gyroradius, SB ~ B, in which case 

., = 2^« 1700 y^(^)"cm. (1) 

In this case the acceleration timescale can be simplified to 

rp:>6(i^)'^^6^, (2) 

where Va is the speed of the shock. For a one-Gauss mag- 
netic field and a relativistic shock {vs c) this gives, for an 
electron with 7 = 10*, acceleration timescale of a few mil- 
liseconds. For lower energy particles the acceleration is even 
faster. With time resolution of the observations of the order 
of minutes, this is "instantaneous" and poses no problems in 
getting the particles to sufficiently high energies. Requiring 
the acceleration rate to be less than variability timescale tvar 
(in the lab frame), requires that 

For protons, however, whose mass and acceleration 
timescale are one thousand times larger than for electrons, 
acceleration can take several minutes. When the observed 
timescales are of the order of minutes, this leads to a situa- 
tion where, for hadronic models, the particle injection spec- 
trum develops during the fiare and has a gradually increas- 
ing maximum energy cutoff even if the observed timescales 
are reduced by strong Doppler boosting. For the remainder 



of this paper we will concentrate on electron acceleration. 
Furthermore, a source of radius 7? cannot confine particles 
with gyro-radius larger than R, so that there is a geomet- 
ric constraint, rg < R, on the maximum particle energy. 
However, this does not place an important constraint on the 
upper-cutoff of the particle distribution for these sources. 

2.2 Converter mechanism 

In the converter process the accelerating particles cross a 
shock from the downstream to the upstream in a neutral 
form, as a neutron or synchrotron photon. This neutral par- 
ticle then decays into a proton and an electron (in the case of 
a neutron) or produces an electron-positron pair (for a pho- 
ton), which can then be scattered again across the shock, 
and the process continues as in the first-order mechanism 
until the accelerating particle enters into the upstream re- 
gion again in a neutral form. In the converter mechanism, 
the average angle at which the upstream particle re-enters 
the shock is larger than for the first-order process, and in the 
converter process a particle can get the maximal F^ boost 
in every cycle instead of just the firs t one, as w as the case 
for the first-order acceleration (Deris hev et al. I [2003; Ster^ 
|200a). 

The scenario involving the shock-crossing in forms of 
free neutrons can probably be excluded from the dominat- 
ing process in this the free neutron lifetime of ~ 15 
minutes, makes the whole cycle too long for these minute- 
scale flares even for F = 50 (corresponding to co-moving 
timescales of a few hours), because multiple shock-crossing 
cycles are needed for the particles to reach sufficient ener- 
gies to account for all of the radiation. For longer duration 
flares, however, even the neutron-based mechanism can be 
plausible. 

Instead, we concentrate on a variant of this process 
where synchrotron photons play the role of the neutral par- 
ticle. This process is expected to work best in ultrarela- 
tivistic shocks where it competes with the first-order mech- 
anism. The converter mechanism is slower at lower ener- 
gies, but soon reaches the rate of the first-order process - 
it can also reach hig her energies than the first-order one 
dPerishev et al.ll2003l ). and thus dominate the acceleration 
of the highest-energy particles. Because both the converter 
mechanism working via synchtrotron photons and the first- 
order process in general work on similar timescales, we ap- 
proximate the converter mechanism to work on timescale 
corresponding to that of the first-order mechanism and in- 
corporate its effect in the first-order electron acceleration 
timescale calculations. 



2.3 Second-order Fermi acceleration 

The second-order, or stochastic, process, accelerates parti- 
cles using scattering centres moving relative to each other 
even without differences in the actual flow speed. For exam- 
ple, Alfven waves in the turbulent downstream of a relativis- 
tic low-Mach-number shock can provid e promising condi- 
tions for efficient stochastic acceleration (|Virtanen fc Vainid 
l2005h with and without a shock. Because the process is not 
tied to the plasma speed, it can continue to accelerate par- 
ticles far away from the shock and for much longer than the 
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first-order process - provided tliere is sufficient turbuience 
present. 

Tfie acceteratio n timescaie for stocfiastic acceleration is 
jRieger et al.ll2007l ) 



Tpil 



VA 



X 



3 cr, 

4 V 



9_ 

2 ' 
A 



wliere tlie Alfven speed, defined by 
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(p + P)/n, witli tlie energy 
^ being a function of tlie 



depends on the enthalpy, h 
density of the plasma, p = nmc 
composition and number density , n. The mass m depends 
on the composition and is mec ~ 2 rric for pure electron- 
positron plasma, and rriep = rric + nip for ionised hydrogen. 
The effect of the gas pressure, P, is taken to be negligible. 

The second order process will be rapid enough to oc- 
cur on a timescaie shorter than the observed flaring time 
provided that 
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where a = 2.1 for an electron-positron plasma and a = 1.9 x 
10'^ for the hydrogen case. In a "highly magnetised plasma" 
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the above constraint for a rapid stochastic process simplifies 
to 



V 104 / V 1 G y V 50 



< 3.5 X 10' 



(8) 



The stochastic acceleration timescaie as a function of 
the magnetic field and the number density of the plasma 
is shown in Fig. [1] for hydrogen plasma and in Fig. [2] for a 
pair plasma. Although the acceleration is very slow when 
the magnetic field is relatively low and the density is high, 
sites such as magnetically dominated AGN jets with rela- 
tively low matter density and compressed magnetic fields 
could favour fast acceleration with the blob matter density 
of the order of ~ 10^-10^ particles (protons and electrons) 
per cm^ and a magnetic field of the order of one gauss, 
thus providing acceleration timescales comparable to the 
observed minute-scale flickering. For purely hadronless pair 
plasma the acceleration timescaie is much shorter, and can 
even be "instantaneous" if the plasma density is not very 
high. In sources where the plasma is purely or mainly lep- 
tonic and has low density, sufficiently high magnetic field 
can turn the second-order acceleration more rapid than the 
first-order one. This, however, requires quite ideal turbu- 
lence conditions with particle-scattering waves moving in 
opposite directions over a sufficiently long lengthscale. The 
Alfven speed is plotted in Fig. |3] for the same plasma pa- 
rameters as in Fig. [T] 

The problem with turbulence conditions in AGN 
sources is that they are very poorly known. This leads 
to large uncertainties regarding the processes that depend 
strongly on the turbulence details - like stochastic accelera- 
tion. In real relativistic sources its efficiency is probably less 
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Figure 1. Stochastic-acceleration timescaie Tpn (Eq. \^ for an 
electron with 7 = 10* as a function of plasma number density n 
and the magnetic field intensity B in ionised hydrogen plasma. 
Timescales vary from years (top left corner) to milliseconds (bot- 
tom right). 
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Figure 2. As Fig.[T] but for pure pair plasma. 



than that given above or simulated by IVirtanen fc Vainiol 
(2005,), because these assumptions exclude the dissipation 
and damping of the turbulence, which would limit the accel- 
eration region spatially. On the other hand, earlier stochas- 
tic acceleration modelling neglects the possibility for much 
stronger turbulence with SB > B as they assume quasilinear 
conditions with 6B <^ B. 



2.4 Shear acceleration 

Fermi acceleration without a shock is also possible where 
ever scattering centres flow at different speeds, even if the 
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Figure 3. Alfven speeds, in the units of c, for parameters corre- 
sponding to those in Fig. [T] 



flows are parallel. An example of this is the longitudinal 
shear across the jet radius, where particles can tap into the 
difference between the fast core of the jet and the slower 
(or even motionless) exterior. This mechanism works both 
on the sharp contact discontinuity between the stream- 
ing jet plasma and the interstellar plasma, as well as on 
smoother flows with gradual shear, with scattering centre 
speed varying slowly as a function of the radius of the jet 
l|Rieger fc Duffvll2004 ). 

The shear acceleration timescale in a relativistic flow, 
with spatial scale J?, depends on the extent of the shear re- 
gion and the scattering timescale according to (|Rieger et al.l 
[20071) 
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With r = 10, B = 1 G, and i? = 10^^ cm we get, for 
an electron with 7 — 10*, timescale of decades and rapid 
acceleration requires that 
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Even for models involvin g minimal size and fast flows 
(see the "needle/jet" model of iGhisellini fc Tavecchicll2008l . 
with i? = 3 X 10^* cm and F = 50), the acceleration is too 
slow to have relevance in minute-timescale flares. Due to 
the strong dependence on the spatial scale (r oc i?^) and the 
jet Lorentz factor (r oc F"*), the shear acceleration requires 
extremely fast and narrow jets or otherwise very sharp tran- 
sition layers between different parts of the jet. 

It is interesting to note, however, that because Tshoar is 
inversely proportional to the particle mean free path and, 
thus, energy - in contrast to the flrst- and second-order 
acceleration - it can, in principle, continue to accelerate 
the highest-energy particles already energised by the flrst- 
or second-order mechanisms. Furthermore, although with 
Bohm-type scattering shear acceleration is too slow, we can 
estimate the circumstances where shear acceleration would 



be signiflcant. We flnd that shear acceleration dominates 
over the first-order process when 
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which is fulfilled for accelerating electrons whose mean free 
path satisfies A > R/F^. For F = 10 this means that the 
mean free path has to be greater than one percent of the size 
of the flaring region. Naturally the mean free path must be 
shorter than the minimum dimension of the blob, leading to 
F~^ < ^ < 1. Depending on the behaviour of the mean free 
path, this condition could be reached easily with different 
scattering models, but as long as we do not have sufficient 
data suggesting such turbulence conditions, we will limit our 
study to Bohm diffusion with A ~ r^. 



3 ENERGY LOSS TIMESCALES 
3.1 Radiation losses 

Previous authors have already addressed the radiation loss 
mechanisms and timescales, so we only review the basic 
limitations the losses put on particle acceleration. A sin- 
gle electron, with Lorentz factor 7, moving in a plasma with 
energy density of the local magnetic field being I7mag and 
the energy density in the radiation field being C/rad, under- 
goes synchrotron and inverse Comptonisation (IC) losses on 
timescales 

where ctt is the Thompson cross section and the energy 
density U in question is I7mag for synchrotron losses and 
[/rad for IC losses. In the latter case. Equation (|12|) assumes 
that Compton scattering happens in the Thomson regime. 
However, this may not be true for high-frequency peaked 
BL Lac sources, such as Markarian 501 and PKS 2155-304 
discussed in the introduction. In the Klein-Nishina regime 
Compton cooling is suppressed by orders of magnitude, and 
the timescale given by Eq. (|12p is radically shorter than what 
would be expected for the highest-energy particles in these 
sources. Detailed cooling time analysis is beyond the scope 
of this acceleration- focused paper, but will be included in 
subsequent studies. 

It is clear from Eq. I|12p that the timescales of the syn- 
chrotron and IC losses differ only in the energy density term. 
For magnetic field the density J/mag = -f^ is straightforward 
to calculate, but the energy density of the radiation field de- 
pends on many things, especially on the source of the pho- 
tons. In this paper we assume a ratio of the energy densities 
and use that as a free parameter. Because of the otherwise 
equal form except for the energy density, the timescales for 
synchrotron and IC losses follow a simple relation: 

Kad 



'synch 
TIC 



(13) 



and IC losses dominate over the synchrotron losses when 
f^rad/f-'^mag > 1, which is Considered to be the case in most 
of the AGN sources with high radiation densities especially 
close to the accretion disc and BLR clouds. Values of the 



order of Kad/t/ir 



lO-" or lO'' are typically used in mod- 



elling these sources. Questions such as whether the dominat- 
ing seed photon field is from an external source or created 
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by the accelerating partic l es (see, e.g. Begelman et al. 20081 : 
ICelotti fc Ghisellinil [2OO8I : iGhisellini fc Tavecchiol I2OO8I . for 

support for both alternatives), or how are the loss timescales 
changed in the extreme Klein-Nishima region, are beyond 
the scope of this paper. Furthermore, we don't include the 
effects of an increas ing magnetic field, but not e that when 
taken into account jSchlickeiser fc Lerche|[2007h the nonlin- 
ear effects due to the change in the equipartition magnetic 
field energy density can lead to synchrotron losses signifi- 
cantly faster than that expected from standard linear losses. 
Even the "nonlinear synchrotron losses" work on timescales 
longer than the fastest acceleration mechanism for the mod- 
est 7 particles. 

3.2 Escape and adiabatic expansion 

For Bohm-type diffusion the electron escape timescale de- 
pends on the size of the acceleration region and the particle 
gyroradius as follows: 
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(14) 



Calculating the timescales it becomes apparent that escape 
can play a significant role in shaping a-few-minute flare only 
if the spatial scale is essentially smaller than what is required 
by the causal connection argument, which limits the emit- 
ting region to sizes smaller than ^ 10^"^ cm (of the order of 
1 AU). Even for spatial scales that small, only particles with 
7 3> fO^ are considerably affected by escape on timescales 
of the order of minutes, so we can safely exclude escaping 
from the dominant restrictions to the particle acceleration 
mechanisms in these short flares. 

We take also the adiabatic losses due to expansion to 
be negligibe within our assumptions. The blob is very likely 
to be travelling within a coUimated flow, where expansion 
is expected to be slow and not to play a major role in the 
first few minutes of the flare. Even though the condition 
and expansion spe ed is not know n for these sources, a rough 
estimate (see, e.g. iLongairllf 9941 ) gives 

-Tad 3 X f —5- s (f5) 

which, in the shortest parsec or sub-parsec scale jets in 
blazars, corresponds to timescales of many hours. For 
sources like microquasars, however, where the speed and the 
degree of coUimation of the jet are probably lower, the adi- 
abatic losses are more likely to set limits for the particle 
energy. 



4 RESULTS & DISCUSSION 

We have compared the various acceleration and loss 
timescales for various different physical environments by 
changing the following parameters: magnetic field intensity, 
matter density, size of the region, the jet Lorentz factor, and 
the ratio of radiation energy density to that of the magnetic 
field. In Fig.|4]the acceleration and loss timescales are shown 
for a s et of parameters that c o rrespo nd to the recent mod- 
els of iGhiseUini fc Tavecchiol (j2008h and iBeeel man et al.l 
l)2008t ) . The bold lines correspond to the fastest acceleration 
(marked with solid lines) and loss (dashed lines) timescales. 
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Figure 4. Calculated timescales in the comoving frame of the 
radiating plasma as a function of the particle Lorentz factor. Ac- 
celeration timescales (solid lines) are shown for the first-order 
Fermi acceleration (labelled Fermi-I) for electron, second-order 
acceleration (Fermi-II) for an electron in a hydrogen plasma, and 
shear acceleration, and losses (dashed lines) due to radiation due 
to synchrotron and IC; escape losses are too slow to appear in 
this Figure. The thick black lines show the fastest gain and loss 
timescales. Values for the magnetic field _B, plasma number den- 
sity n, bulk Lorentz factor F, and the ratio J7rad/'^mag are given 
in the figure. 



and their crossing point corresponds to the Lorentz factor 
where gains and losses are equal - beyond this particles cool 
instead of accelerating. The most severe energy loss mecha- 
nism is the inverse Compton process, due to L^rad/Knag > 1. 
With the aforementioned assumptions, the particle escape 
process is orders of magnitude too slow to affect the radi- 
ating electrons significantly within timescales of minutes or 
hours. 

The spatial scale affects the escape and shear accelera- 
tion timescales in the same way: both processes become more 
important for smaller sizes, and when considering minute- 
or even hour-scale variability in sources with R ~ 10^^ cm, 
both are negligible. It is possible, however, that if the parti- 
cle scattering mean free path is significantly larger than the 
particle gyroradius, both processes begin to have noticeable 
effects. 

The first-order process (including the effects of con- 
verter mechanism) accelerates electrons up to 7 ~ 10*' in 
just seconds, so it is fair to say that even in these short 
flares acceleration can be considered instantaneous. How- 
ever, at the same time the second-order mechanism works 
on slower but still interesting timescales. 

The plasma density directly affects only the second- 
order acceleration rate through the Alfven-speed depen- 
dence: higher density leading to slower acceleration. Compo- 
sition of the plasma has the same effect, because the presence 
of heavier ions decreases the Alfven speed. 

The expected rapid second-order acceleration process 
raises two important questions. Firstly, is the effect of 
stochastic acceleration seen in the observations? And, sec- 
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Figure 5. An example of parameter set which leads to a possi- 
ble dominance of acceleration at all energies: low magnetic field 
and strong shear make the shear acceleration energise high-energy 
particles faster than the IC radiation decreases the energy. 

ondly, if it is not an important process, how do we ex- 
plain this? For the first question, some recent observations 
of particle spectra with very hard power-law spectral in- 
dices (a < 2, for A'"(7 ) oc 7"°^) suRgest w e may see it 
iKatarzvnski et alJlioog : iBottcher et alj|2007l '). 

The latter question, i.e., if we do not see stochastic ac- 
celeration even if the conditions seem perfect for it, suggests 
basically three plausible options. 

Turning to the second question, an absence of stochastic 
acceleration would raise the following possibilities. Firstly, 
the acceleration model could be wrong. However, for the rea- 
sons given later in the discussion, we do not believe this to be 
the case and leave this option for other authors. Secondly, it 
is likely that we do not fully understand the plasma param- 
eters such as composition and the magnetic field. Thirdly, 
the conditions are right, but there are other factors suppress- 
ing the acceleration or turning it off. This is also a realis- 
tic scenario, as we have very little knowledge of the nature 
and extent of the turbulence in the acceleration and radiat- 
ing regions. A turbulence damping process could shut down 
this mechanism. This option is very plausible, because the 
process depends heavily on the turbulence in the accelera- 
tion regions, and very little is known about it. Furthermore, 
the stochastic acceleration in relativistic shocks and AGN- 
related conditions is still little studied. It is not unlikely that 
future studies will be able to set stricter limits for the ex- 
tent (both spatial and energetic) of the process, making only 
a small region immediately behind the turbulence-creating 
shock front suitable for stochastic acceleration. 

4.1 Absence of continuous acceleration 

It is interesting to note that with the assumption of high- 
velocity plasma streaming (whether as a continuous jet or 
an ejected plasma "blob") within an external plasma, there 
always seems to be suitable conditions for very rapid particle 



acceleration. The acceleration can be "instantaneous" first- 
order acceleration - feeding the radiating region constantly 
with a power law of high-energy particles - or it can be a 
more gradual process, changing the energy distribution of 
the radiating particles on longer timescales. 

When looking for a set of parameters that would yield 
dominating energy losses - corresponding to models with in- 
stantaneous acceleration to a power-law spectrum followed 
by only radiative cooling and other losses - we could not 
find any within limits typically considered physical for these 
sources: B € [W~^,l] G; F £ [2,50]; n G [1,10^°] cm'^ 
R G [10^", lO^'^] cm; Kad/C/mag G [10-^10=*]. Acceleration 
times within these limits were always shorter than the loss 
timescales in the lower energies up to electron Lorentz fac- 
tors 7 ~ 10^~^. In fact, in some cases the "inverse" energy 
dependence of the shear acceleration led to domination of 
acceleration on all energies - an example of this is shown in 
Fig-E] Although we don't expect this to be a typical case in 
real sources (for reasons discussed earlier) it underlines the 
problem of getting only cooling without having acceleration 
at the same time. 

4.2 Time delays 

Interesting constraints on acceleration mechanisms arise 
from observations of a TeV fiare of Markarian 501 on July 
9th, 2005. Here the hard gamma-rays (2.7 TeV) were ob- 
served to lag behind the soft ones (190 GeV) by roughly 
four minutes. Similar "hard lags " have been observed ea rlier 
also in the X-ray domain (e.g., iBrinkmann et aL I l2005l . for 
Markarian 401), but no definite explanation has been found 
yet. The "soft lags" - where lower-energy radiation peaks 
later than the higher-energy one - are easier to explain by 
high-energy radiating particle cooling and radiating on lower 
and lower frequencies, but, intuitively, these hard lags would 
require particles not cooled but heated or accelerated during 
the fiare. 

The observed delay has been proposed as being sue 
to the rad iating blob accelerating during the flare by 
iBednarek fc Wagner (2008), but even in their model the 
particles would only undergo cooling without any ongo- 
ing acceleration in or around the high-F plasma flow. 
iMastichiadis fc MoraitisI l|2008l ), however, showed that the 
observed features can be explained with a very simple and 
physically plausible modification, namely allowing the par- 
ticles to accelerate gradually. In their model the acceleration 
time is of the order of hours and the electrons reach energies 
corresponding to 7 ~ 10®. The timescale may be too long for 
the first-order mechanism, but, interestingly, matches well to 
the stochastic acceleration timescale corresponding to their 
parameters. 

Furthermore, similar slow en ergisation was quantita- 
tively predicted bv lVirtanen fc Vai nio (2005), whose simula- 
tions of stochastic acceleration in relativistic shocks showed 
a gradual shift of the whole particle spectrum to higher en- 
ergies as the particles were accelerated behind a shock front. 
Hadronic process could also be relevant within this model. 

4.3 Generation of turbulence 

Finally, one must also remember that as the accelerating 
particles need electromagnetic turbulence for scattering and 
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energy gain, it is obvious that suitable turbulence either has 
to exist in the pre-shock plasma or be generated at the shock 
by the particles themselves. All models with instantaneous 
or constant first-order Fermi acceleration require the acceler- 
ation timescale to be significantly less than the flaring time, 
and the same limit is set also for the turbulence generation. 
Recent numerical work has addressed the problem of turbu- 
lence generation, and although the pr oblem is still far from 
being solved, studies such as that by iReville et all (|2006t ). 
suggest that especially nonlinear t urbulence generation by 
the particles themselves (found bv iBelll [20041 ) could indeed 
be efficient enough to provide sufficient turbulence for fast 
acceleration. 



5 CONCLUSIONS 

We have estimated the timescales of different particle- 
acceleration mechanisms in the context of minute-scale TeV 
blazar flares. We have excluded both the neutron-based con- 
verter mechanism and shear acceleration from the dominat- 
ing processes in these sources on observed timescales of the 
order of minutes. Furthermore, the flrst-order Fermi accel- 
eration of protons is likely to be too slow to appear "in- 
stantaneous". Instead, the accelerated proton distribution 
is expected to continue to be energised during the flare and 
could have special signiflcance in hard-lag sources. 

Various simultaneously active acceleration mechanisms 
working on different timescales suggest that in these sources 
one could expect to find slower, gradual energisation in ad- 
dition to instantaneous - or a brief initial period of - accel- 
eration. This speaks in favour of models with ongoing accel- 
eration, instead of models with instantaneous injection of a 
fully-developed power-law spectrum which then only under- 
goes cooling. However, since many models that only include 
coofing can reproduce the observations, their results com- 
bined with timescale analysis could be helpful in studying 
the source parameters. 

We emphasise that the current model is very simpli- 
fied. The results neglect, in particular, the Klein-Nishina 
effects in the energy losses of the highest-7 particles. Also, 
the magnetic field is assumed to be constant throughout the 
fiare, although generation of turbulence and the subsequent 
increase in the magnetic field strength can decrease the syn- 
chrot ron loss timescales significantly (|Schlickeiser fc Lerchd 
I2OO7I ). Furthermore, because the acceleration efficiency de- 
pends on the scattering mean free path, for different scat- 
tering models the acceleration timescales could also change. 
Furthermore, other omitted turbulence-affecting effects can 
decrease the acceleration timescales, for example by increas- 
ing the acceleration efficiency. An example is the case of 
turbulence transmission leading to an increased scattering- 
centre compressio n rati o (see Vainio et al., 2003; 20051 . and 
Tammi & Vainio 120061 ) ■ and enhanced first-order acceler- 
ation at th e shock in co nditions similar to AGN or micro- 
quasar jets (|Tammill2008l ). However, even in its current state, 
especially in the lower-energy part of the particle spectrum 
(7 ^ 10*) the present analysis still provides a useful tool for 
estimating the source properties. 

Finally, we note that since the acceleration efficiency 
depends on the scattering mean free path differently for each 
of the processes, detailed modelling could provide additional 



tests when we learn more about the intrinsic properties of 
the plasma and turbulence in these sources. 
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